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In this article, we investigate the contribution of the higher- twist Feynman diagrams to the 
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coupling approaches. The structure of infrared renormalon singularities of the higher twist subpro- 
cess cross section and the resummed expression (the Borel sum) for it are found. We compared the 
resummed higher-twist cross sections with the ones obtained in the framework of the frozen coupling 
approach and leading-twist cross section. We obtain, that ratio R = (E^+) res /(^^+)°, for all val- 
ues of the transverse momentum px of the pion identically equivalent to ratio r=(A^ T ) res / (A^ T )°. 
It is shown that the resummed result depends on the choice of the meson wave functions used in 
calculation. Phenomenological effects of the obtained results are discussed. 
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I. INTRODUCTION 



The large-order behavior of a perturbative expansion in gauge theories is inevitably dom- 
inated by the factorial growth of renormalon diagrams [1-4]. In the case of quantum chro- 
modynamics (QCD), the coefficients of perturbative expansions in the QCD coupling a s can 
increase dramatically even at low orders. This fact, together with the apparent freedom in 
the choice of renormalization scheme and renormalization scales, limits the predictive power 
of perturbative calculations, even in applications involving large momentum transfers, where 
a s is effectively small. 

A number of theoretical approaches have been developed to reorganize the perturbative 
expansions in a effort to improve the predictability of the perturbative QCD (pQCD). For 
example, optimized scale and scheme choices have been proposed, such as the method of 
effective charges (ECH) [5], the principle of minimal sensitivity (PMS) [6], and the Brodsky- 
Lepage- Mackenize (BLM) scale-setting prescription [7] and its generalizations [8-20]. In [4] 
developments include the resummation of the formally divergent renormalon series and the 
parametrization of related higher-twist power-suppressed contributions. 

In general, a factorially divergent renormalon series arises when one integrates over the 
logarithmically running coupling a s (k 2 ) in a loop diagram. Such contributions do not occur 
in conformally invariant theories which have a constant coupling. Of course, in the physical 
theory, the QCD coupling does run. 

Among the fundamental predictions of QCD are asymptotic scaling laws for large-angle 
exclusive processes [21-27]. QCD counting rules were formalized in Refs. [22,23]. These re- 
actions probe hadronic constituents at large relative momenta, or equivalently, the hadronic 
wave function at short distances. Important examples of exclusive amplitudes are provided 
by the electromagnetic form factors of mesons. Since there is little direct evidence with 
which to compare the predictions, it is fortunate that short-distance wave functions also 
control a wide variety of processes at large transverse momentum. In particular, the meson 
wave function determines the leading higher-twist contribution to meson production at high 
Vt- 

The hadronic wave functions in terms of quark and gluon degrees of freedom play an im- 
portant role in the quantum chromodynamics predictions for hadronic processes. In the per- 
turbative QCD theory, the hadronic distribution amplitudes and structure functions which 
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enter exclusive and inclusive processes via the factorization theorems at high momentum 
transfers can be determined by the hadronic wave functions, and therefore they are the 
underlying links between hadronic phenomena in QCD at large distances (nonperturbative) 
and small distances (perturbative). If the hadronic wave functions were accurately known, 
then we could calculate the hadronic distribution amplitude and structure functions for 
exclusive and inclusive processes in QCD. Conversely, these processes also can provide phe- 
nomenological constraints on the hadronic distribution amplitudes, the hadronic structure 
functions, and thereby the hadronic wave functions. From another point of view, as different 
wave functions may give the same distribution amplitude, there are still ambiguities about 
the wave function even if we know the exact form of the distribution amplitude. 

In principle the hadronic wave functions determine all properties of hadrons. From the 
relation between the wave function and measurable quantities we can get some constraints 
on the general properties of the hadronic wave functions. Note that only the lowest qq 
Fock state contributes to the leading scaling behavior; other Fock-state contributions are 
suppressed by powers of l/Q 2 . 

The frozen coupling constant approach can be applied for investigation, not only exclu- 
sive processes, but also for the calculation of higher-twist contributions to some inclusive 
processes, for example as large -pr meson photoproduction [28], two-jet +meson production 
in the electron-positron annihilation [29]. In the works [28,29] for calculation of integrals, 
such as 



the frozen coupling constant approach was used. According to Ref . [7] should be noted that in 
pQCD calculations the argument of the running coupling constant in both, the renormaliza- 
tion and factorization scale Q 2 should be taken equal to the square of the momentum transfer 



x\ — > 0, x 2 — > integrals (1.1) diverge and for their calculation some regularization methods 
of a s (Q 2 ) in these regions are needed. The power-suppressed corrections arising from the 
soft end-point regions to the single meson photoproduction process were computed in[31]. 
In Ref. [31], the evolution of the meson wave function on the factorization scale was ignored. 
In the present work, we take into account this evolution as well. In Ref. [32], the authors in- 




(1.1) 



of a hard gluon in a corresponding Feynman diagram. But defined in this way, a s (Q 2 ) suf- 
fers from infrared singularities. For example in our work [30], Q 2 equals to (xi — l)u and 
—Xit, where u, i are the subprocess's Mandelstam invariants. Therefore, in the soft regions 
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vestigated the phenomenology of infrared renormalons in inclusive processes. The dispersive 
approach has been devised to extend properly modified perturbation theory calculations 
towards the low-energy region [33]. Connections between power corrections for the three 
Deep Inelastic Scattering sum rules have also been explored in [34]. 

Investigation of the infrared renormalon effects in various inclusive and exclusive pro- 
cesses is one of the most important and interesting problems in the perturbative QCD. 
It is known that infrared renormalons are responsible for factorial growth of coefficients 
in perturbative series for the physical quantities. But, these divergent series can be re- 
summed by means of the Borel transformation [1] and the principal value prescription [35], 
and effects of infrared renormalons can be taken into account by a scale-setting procedure 
a s{Q 2 ) — ¥ a s{. ex v(f (Q 2 ))Q 2 ) & t the one-loop order results. Technically, all-order resum- 
mation of infrared renormalons corresponds to the calculation of the one-loop Feynman 
diagrams with the running coupling constant a s (—k 2 ) at the vertices or, alternatively, to 
calculation of the same diagrams with nonzero gluon mass. Studies of infrared renormalon 
problems have also opened new prospects for evaluation of power-suppressed corrections 
to processes characteristics [36]. Power corrections can also be obtained by means of the 
Landau-pole free expression for the QCD coupling constant. The most simple and elabo- 
rated variant of the dispersive approach, the Shirkov and Solovtsov analytic perturbation 
theory, was formulated in Ref.[37]. 

By taking these points into account, it may be argued that the analysis of the higher-twist 
effects on the dependence of the pion wave function in pion production at proton-proton 
collisions by the running coupling (RC) approach [38], are significant from both theoretical 
and experimental points of view. 

We will show that higher-twist terms contribute substantially to the inclusive meson cross 
section at moderate transverse momenta. In addition, we shall demonstrate that higher- 
twist reactions necessarily dominate in the kinematic limit where the transverse momentum 
approaches the phase-spase boundary. 

A precise measurement of the inclusive charged pion production cross section at \/s = 
62.4 GeV and y^s = 200 GeV is important for the proton-proton collisions program at 
the Relativistic Heavy Ion Collider (RHIC) at the Brookhaven National Laboratory. The 
result of our calculations are based on the proton-proton collisions at y^s = 62.4 GeV and 
y/s = 200 GeV 
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Another important aspect of this study is the choice of the meson model wave functions. 
In this respect, the contribution of the higher-twist Feynman diagrams to a pion production 
cross section in proton-proton collisions has been computed by using various pion wave 
functions. Also, higher-twist contributions which are calculated by the running coupling 
constant and frozen coupling constant approaches have been estimated and compared to 
each other. Within this context, this paper is organized as follows: In SecJTTl we provide 
formulas for the calculation of the contribution of the high twist diagrams. In Sec. II I II we 
present formulas and an analysis of the higher-twist effects on the dependence of the pion 
wave function by the running coupling constant approach. In Sec. IIV} we provide formulas 
for the calculation of the contribution of the leading-twist diagrams. In Sec. EJ we present 
the numerical results for the cross section and discuss the dependence of the cross section 
on the pion wave functions. We state our conclusions in Sec. IVII 

II. CONTRIBUTION OF THE HIGH TWIST DIAGRAMS 

The higher-twist Feynman diagrams, which describe the subprocess q\ + q 2 — > 7r + (7r~) +7 
for the pion production in the proton-proton collision are shown in Fig.l. In the higher-twist 
diagrams, the pion of a proton quark is directly observed. Their 1/Q 2 power suppression 
is caused by a hard gluon exchange between pion constituents. The amplitude for this 
subprocess can be found by means of the Brodsky-Lepage formula [26] 

M(s,t)= / dx 1 / dx 2 5(l -x 1 - x 2 )§ 7T (x 1 ,x 2 ,Q 2 )T H (s, t;x 1 ,x 2 ). (2.1) 
Jo Jo 

In Eq.(2.1), Tjj is the sum of the graphs contributing to the hard-scattering part of the 
subprocess. The hard-scattering part for the subprocess under consideration is q\ + q 2 — > 
(9192) + 7, where a quark and antiquark form a pseudoscalar, color-singlet state (qiq 2 ). 
Here §(xi,x 2 ,Q 2 ) is the pion wave function, i.e., the probability amplitude for finding the 
valence qiq 2 Fock state in the meson carry fractions X\ and x 2 , X\ + x 2 = 1. Remarkably, 
this factorization is gauge invariant and only requires that the momentum transfers in Tjj 
be large compared to the intrinsic mass scales of QCD. Since the distribution amplitude and 
the hard-scattering amplitude are defined without reference to the perturbation theory, the 
factorization is valid to leading order in 1/Q, independent of the convergence of perturbative 
expansions. 
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The hard-scattering amplitude T# can be calculated in perturbation theory and rep- 
resented as a series in the QCD running coupling constant a s (Q 2 ). The function $ is 
intrinsically nonperturbative, but its evolution can be calculated perturbatively. 

The q{q 2 spin state used in computing T H may be written in the form 



where e± = =f(1/ v / 2)(0, 1, ±i, 0) in a frame with (pm)i,2 = and the N* project out 
a state of a spin s, and pu is the four-momentum of the final meson. In our calculation, 
we have neglected the pion and proton masses. Turning to extracting the contributions 
of the higher-twist subprocesses, there are many kinds of leading-twist subprocesses in pp 
collisions as the background of the higher-twist subprocess q\ + q 2 — > n + (or 7r~) +7, such 
as q + q -> 7 + g(g -> 7r + (7r-)), g + g -> 7 + g(g -> 7r + (7r-)), g + 5 -> 7 + g#(g -> 
7r + (7r - )) etc. The contributions from these leading- twist subprocesses strongly depend on 
some phenomenological factors, for example, quark and gluon distribution functions in the 
proton and fragmentation functions of various constituents, etc. Most of these factors have 
not been well determined, neither theoretically nor experimentally. Thus they cause very 
large uncertainty in the computation of the cross section of process pp — > ix + (or 7T~) +7 + X. 
In general, the magnitude of this uncertainty is much larger than the sum of all the higher- 
twist contributions, so it is very difficult to extract the higher-twist contributions. 

The Mandelstam invariant variables for subprocesses q 1 + q 2 — > it + (it~) + 7 are defined 

as 



In our calculation, we have also neglected the quark masses. We have aimed to calculate 
the pion production cross section and to fix the differences due to the use of various pion 
model functions. We have used five different wave functions: the asymptotic wave function 
(ASY), the Chernyak-Zhitnitsky wave function [27,39], the wave function in which two 
nontrivial Gegenbauer coefficients a 2 and 04 have been extracted from the CLEO data on 
the 77* — > 7T° transition form factor [40], the Braun-Filyanov pion wave functions [41] and 
the Bakulev-Mikhailov- Stefanis pion wave function[42]. The wave functions of pions also are 
developed in Refs. [43-45] by the Dubna group. In Ref.[40], the authors have used the QCD 




(2.2) 




(2.3) 
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light-cone sum rules approach and have included into their analysis the NLO perturbative 
and twist-four corrections. They found that in the model with two nonasymptotic terms, at 
the scale fx = 2.4 GeV, a 2 = 0.19, a 4 = -0.14. 

$asy(x) = VSf w x(l - x), <f> C z(x, l4) = ^asy^X - l) 2 , 

®cleo(x, nl) = $a Sy (x)[l + 0.405(5(2x - l) 2 - 1) - 0.4125((2x - l) 4 - U(2x - l) 2 + 1), 

$ BF (x, nl) = ®asy(x)[l + 0.66(5(2x - l) 2 - 1) + 0.4687((2:r - l) 4 - U(2x - l) 2 + 1)], 

<5>bms{x, f4) = W*) [1 + 0.282(5(2x - l) 2 - 1) - 0.244((2:r - l) 4 - 14(2x - l) 2 + 1)], (2.4) 

where f n = 0.923GeV is the pion decay constant. Here, we have denoted by x = xi, 
the longitudinal fractional momentum carried by the quark within the meson. Then, x 2 = 
1 — x and xi — x 2 — 2x — 1. The pion wave function is symmetric under the replacement 
xi — x 2 <-> x 2 — x\. The model functions can be written as 



$asy(x) = V3f v x(l -X), 



®Cz(x,fl 2 Q ) = <5> asy (x) 



C 3 /2 (2x-l) + ~-Cr(2x-l) 

,3/2 



,3/2, 



CLEO 



[x^D 



$asy(x) 



Cl /2 (2x - 1) + 0.27C2 /Z (2a; - 1) - $.22Ct'\2x - 1) 



3/2, 



$ BF (x,/i 2 ) = <S> asy (x) \cl ,2 (2x - 1) + 0.44C 2 3/2 (2:r - 1) + 0.25C 4 3/2 (2:r - 1) 



$ bms(x, V 2 ) = <$> asy (x) Cl /2 (2x - 1) + 0.188C* /Z (2a; - 1) - 0.13Cf \2x - 1) 



3/2, 



,3/2, 



Cl'\2x - 1) = 1, C*<\2x - 1) = ^(5(2* - l) 2 - 1), 



,3/2, 



Cl /2 (2x - 1) = ^(21(2x - l) 4 - 14(2x - l) 2 + 1). 
8 

Several important nonperturbative tools have been developed which allow specific pre- 
dictions for the hadronic wave functions directly from theory and experiments. The QCD 
sum-rule technique and lattice gauge theory provide constraints on the moments of the 
hadronic distribution amplitude. As is seen from Eq.(2.5) wave functions of meson, which 
were constructed from theory and experiment strongly depend on the methods is applied. 
However, the correct pion wave function is still an open problem in QCD. It is known that the 
pion wave function can be expanded over the eigenfunctions of the one-loop Brodsky-Lepage 
equation, i.e., in terms of the Gegenbauer polynomials {Cn 2 {2x — 1)}, 
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(2.5) 
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$ v (x,Q 2 ) = <S> asy {x) 



1+ £ *n{Q 2 )Cl /2 {2x-l) 



n=2,4.. 



(2.6) 



The evolution of the wave function on the factorization scale Q 2 is governed by the 
functions a n (Q 2 ), 



a n (Q ) = On(^o) 



a«(A*o) 



(2.7) 



72 = 50 74 = 364 
A) 81' A> 405' % 
In Eq.(2.7), {7„} are anomalous dimensions defined by the expression, 



In = Ci 



n+1 n 



(n + l)(n + 2) 



(2.8) 



The constants a n (//o) = a n are input parameters that form the shape of the wave functions 
and which can be extracted from experimental data or obtained from the nonperturbative 
QCD computations at the normalization point //q. The QCD coupling constant a s (Q 2 ) at 
the one-loop approximation is given by the expression 



a s (Q 2 ) 



An 



(2.9) 



A,/n(Q 2 /A 2 )- 

Here, A is the fundamental QCD scale parameter, /3q is the QCD beta function one-loop 
coefficient, respectively, 

2 

A) = ll - gW/. 

The cross section for the higher-twist subprocess q\q 2 — > is given by the expression 

r 1 1 1 

(2.10) 



da,, f ^ 8n 2 a E C F [D(t,u)\ 
^ (a '*' tt) = — 27 ¥~ 



u 



1 1 

+ ¥ 



where 



D(i, u) = e\i l dx\ 
Jo 



a s (Ql)®n(xi,Ql) 



+ 6 2 U I dX\ 

'o 



1 — X\ 



(2.11) 



Here Q\ = (xi — l)u, and Q\ = —x±i, represent the momentum squared carried by the 
hard gluon in Fig.l, ei(e 2 ) is the charge of qi(q 2 ) and Cf = f ■ The higher-twist contribution 
to the large- px pion production cross section in the process pp — > + 7 is [46] 



d(j f f s da 

E m T = E-^-= / / dx 1 dx 2 G qi/hl (x 1 )G q2/h2 (x 2 )-^(qq^7r-f)5(s + i + u). (2.12) 
« V Jo Jo 71 dt 



^ da da 
irE 



d 3 p dydp\ ' 

S = X±X 2 S, 
t = Xit, 

u = x 2 u, (2-13) 



t = —rriTy/se y = —pr^/se y , 

u = —m T \fse y = —p T y/se y , 
x 2 u x 2 p Ty /se y 



x 1 



x 2 



x 2 s + t x 2 s — pTy/se y ' 

X\t X\PT\J~SZ~ y 



x\s + u x\s — PT\fse y '' 
where my - is the transverse mass of pion, which is given by 

2 2,2 
m T = m + p T . 

Let us first consider the frozen coupling approach. In this approach we take the four- 
momentum square Q\ 2 of the hard gluon to be equal the pion's transverse momentum square 
Q\,2 = Pt- m ^is case, the QCD coupling constant a s in the integral (2.11) does not depend 
on the integration variable. After this substitution calculation of integral (2.11) becomes 
easy. Hence, the effective cross section obtained after substitution of the integral (2.11) into 
the expression (2.10) is referred as the frozen coupling effective cross section. We will denote 
the higher-twist cross section obtained using the frozen coupling constant approximation by 

For a full discussion, we consider a difference A^ T between the higher-twist cross section 
combinations H H T and Y, HT 

A HT = sf T _ sf _T = E^{pp - 7T+7) - E^{pp - ,r- 7 ). (2.14) 

We have extracted the following higher-twist subprocesses contributing to the two covari- 
ant cross sections in Eq.(2.12) 

^(wd^TT+7), ^(Jm^TT+7), ^(Mrf^7T"7), ^(rfM^7T^7), (2.15) 

By charge conjugation invariance, we have 

da 1 - , da 3 _ da 2 - , da 4 _ 

— —{ud — > 7r 7) = — — (ud — > 7r 7), and — —(du — > ir 7) = — —(du — > n 7). (2.16) 
dt dt dt dt 
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III. THE RUNNING COUPLING APPROACH AND HIGHER-TWIST MECHA- 
NISM 

In this section we shall calculate the integral (2.11) using the running coupling constant 
method and also discuss the problem of normalization of the higher-twist process cross 
section in the context of the same approach. 

As is seen from (2.11), in general, one has to take into account not only the dependence 
of ot{Q\ 2 ) 011 the scale Q\ 2 , but also an evolution of $(x, Q\ 2 ) with Q\ 2 . The meson wave 
function evolves in accordance with a Bethe-Salpeter-type equation. Therefore, it is worth 
noting that, the renormalization scale (argument of a s ) should be equal to Q\ = (xi — l)u, 
Q\ = —xit, whereas the factorization scale [Q 2 in $m(^, Q 2 )] is taken independent from x, 
we take Q 2 = "p\. Such approximation does not considerably change the numerical results, 
but the phenomenon considered in this article (effect of infrared renormalons) becomes 
transparent. The main problem in our investigation is the calculation of the integral in (2.11) 
by the running coupling constant approach. The integral in Eq.(2.11) in the framework of 
the running coupling approach takes the form 



The a s (XQ 2 ) has the infrared singularity at x — > 1, if A = 1 — x and as a result integral 
(3.1) diverges (the pole associated with the denominator of the integrand is fictitious, because 
$m ~ (1 — x), and therefore, the singularity of the integrand at x — 1 is caused only by 
a s ((l — x)Q 2 )). For the regularization of the integral we express the running coupling at 
scaling variable a s (XQ 2 ) with the aid of the renormalization group equation in terms of the 
fixed one a s (Q 2 ). The renormalization group equation for the running coupling a = a s /ir 
has the form [35] 




(3.1) 



da(XQ 2 ) 



f[«(AQ 2 )] 2 



(3.2) 



dlnX 



where 



A) = 11 - gW/. 



The solution of Eq.(3.2), with the initial condition 



a(X)\ x =i = a = a s (Q 2 )/n, 
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is [35] 



a(X) 



a 



1 + a — InX 
4 



(3.3) 



This transcendental equation can be solved iteratively by keeping the leading a k ln k \ 
order. This term is given by 

a(XQ 2 ) 



1 + InX/t 

After substituting Eq.(3.4) into Eq.(2.11) we get 



(3.4) 



D(Q 2 ) = eit / dx 



a s (XQ 2 )<S> M (x,Q 2 ) 



1-x 



i 

+ €2U I dX\ 




a s {XQ 2 )<b M {x,Q 2 ) 



■ fn^ l\ ®m(x,Q 2 ) a 2 Z" 1 

eiia s (C^ ) / dxi— — — — + e 2 ua s {Q ) / da; 



o 1 (1 - x)(l + InX/t) 



1-x 
<5>m{x,Q 2 ) 



o (l-x)(l + ln\/t) 



e\ta s {Q ) / <ir- 
Jo 



1 $asy(x) 



1 + 5^2,4,.. a n(/^0 



«.(Q') 



7n/0O 



Cl /2 (2x - 1) 



(1 -x)(l +ln\/t) 



e\ua s (Q ) / <ix- 
Jo 

where t = 47r/a: s ((5 2 )A) 



1 $asy(x) 



1 + X^2,4,.. a n(/^C 



q.(Q') 



Cl ,2 {2x-l) 



[1 - x)(l + InX/t) 



(3.5) 



The integral (3.5) is common and, of course, still divergent, but now it is recast into a 
form, which is suitable for calculation. Using the running coupling constant approach, this 
integral may be found as a perturbative series in a s (Q 2 ) 



n=l ^ ' 



(3.6) 



The expression coefficients S n can be written as power series in the number of light quark 
flavors or, equivalently, as a series in power of (3q. 



n-l 



S n — C u /3q 

The coefficients C n of this series demonstrate factorial growth C n ~ (n — 1)!, which might 
indicate an infrared renormalon nature of divergences in the integral (3.5) and corresponding 
series (3.6). The procedure for dealing with such ill-defined series is well known; one has to 
perform the Borel transform of the series [15] 



B[D](u) = J2 



-u 



n=0 
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then invert B[D](u) to obtain the resummed expression (the Borel sum) D(Q 2 ). After this 
we can find directly the resummed expression for D(Q 2 ). The change of the variable x to 
z = ln(l — x), as ln{l — x) = InX. Then, 



D{Q 2 ) = ei ta s (Q 2 )t 



1 ® M (x,Q 2 )dx 



+ e 2 ua s (Q 2 )t 



1 $ M (x,Q 2 )dx 
(l-x)(t + z) 



(3.7) 



/o (l-x)(t + z) 

For the calculation of the expression (3.7) we will apply the inverse Laplase transform to 
Eq.(3.7) [47]. After this operation, formula (3.7) is simplified and we can extract the Borel 
sum of the perturbative series (3.6) and the corresponding Borel transform in dependence 
from the wave functions of the meson, respectively. Also after such manipulations the 
obtained expression can be used for numerical computations. 

The inverse Laplace transformation from l/(t + z) has the form: 



— =r 

t + z J 



e-^ t+z)u du 



(3.8) 



after inserting Eq.(3.8) into (3.7). Then, we obtain 



D(Q 2 ) = ei ia s (Q 2 )t 
e 2 ua s (Q 2 )t 



1 poo 

Jo 



<& M (x,Q 2 )e~ {t+z)u dudx 



1 poo 

Jo 



(1-x) 

§ M (x,Q 2 )e-( t+z ^ u dudx 
O 3 ^ ' 



(3.9) 



In the case of & asy (x) for D(Q 2 ), we get 
D(Q 2 ) - 



4v / 37r/ 7r ei£ 4^/3n f w e 2 u 



A) A) 

In the case of the &cz(x, Q 2 ) wave function, we find 



due 



-tu 



1 — u 2 — u 



(3.10) 



D{Q 2 ) 



4v / 37r/ 7r eit 4v / 37r/ 7r e 2 M 



A) 



due 



—tu 



0.84 



MQ 2 y 


50/81 r 


m a a (t4). 





24 40 
+ 



1 — u 2 — u 
20 



1 — u 2 — u 3 — u 4 — u 
In the case of the <&cleo(x, Q 2 ) wave function, we get 



D (q2^ = f 4y / 37r/ 7r e 1 t + 4 v / 37r/ 7r e 2 M 



A) 



A 



24 40 
+ 



due 



20 



- tu 



1 



1 



1 — u 2 — u 



+ 0.405 



g s (Q 2 ) 

a s {n 2 ) 



(3.11) 



50/81 



1 — u 2 — u 3 — u 4 — if 



- 0.4125 



« s (g 2 ) 



364/405 



12 



120 560 

+ 



1112 1008 

+ 



336 



1 — u 2 — u 3 — u 4 — u 5 — u 6 — u 
Also, in the case of the Q 2 ) wave function, we get 



D{Q 2 ) = 



4v / 37r/ 7r eit 4v / 37r/ 7r e 2 u 



A) 



A) 



24 40 
+ 



Jo 



due 



20 



—in 



1 



1—u 2—u 



■f 0.282 



a s (Q 2 ) 



(3.12) 



50/81 



1 — u 2 — u 3 — u 4 — u 



- 0.244 



120 560 

+ 



1112 1008 

+ 



a s (Q 2 ) 
336 



364/405 



(3.13) 

1 — u 2 — u 3 — u 4 — u 5 — u 6 — u 
Equation(3.1) and (3.2) is nothing more than the Borel sum of the perturbative series (3.6), 
and the corresponding Borel transform in the case $ asy (x) is 



(3.14) 



in the case <&cz(x, Q 2 ) is 



1 



B[D](u) = — 



1 



1-0.84 



u 2 — u 
in the case &cleo(x, Q 2 ) is 



a s (Q 2 



2\\ 50/81 



24 40 
+ 



20 



1 — u 2 — u 3 — u 4 — u 



, (3.15) 



B[D](u) = 



1 



0.4125 



1—u 2—u 

/ n2 v\ 364/405 



<Xs((J>o) 



+- 0.405 

8 



g s (Q 2 

« s («o) 
120 



2\\ 50/81 



24 40 
+ 



20 



+ 



560 



1 — u 2 — u 3 — u 4 — u 

1112 1008 336 
+ 



1 — u 2 — u 3 — u 4 — u 5 — u 6 — u 



(3.16) 



and in the case <&bms(x, Q 2 ) is 
1 1 



B[D\{u) 



0.244 



1—u 2—u 
g s (Q 2 ^ 364/405 



+ 0.282 



g s (Q 2 

«s(«o) 



2 \\ 50/81 



24 40 
+ 



20 



120 560 

+ 



1 — u 2 — u 3 — u 4 — u 

1112 1008 336 
+ 



1 — u 2 — u 3 — u 4 — u 5 — u 6 — u 



(3.17) 



The series (3.6) can be recovered by means of the following formula 

(j \ n— 1 
Tu) fl ™«>l~ 

The Borel transform B[D\{u) has poles on the real u axis at u = 1; 2; 3; 4; 5; 6, which confirms 
our conclusion concerning the infrared renormalon nature of divergences in (3.6). To remove 
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them from Eqs.(3.11) and (3.12) some regularization methods have to be applied. In this 
article we adopt the principal value prescription. We obtain: in the case $ as „ 



[D(Q 2 )r 

in the case & C z(x, Q 2 ) 
[D(Q 2 )} res = 



4 v / 37r/ 7r ei£ 4 v / 37r/ 7r e 2 M 



A 



4v / 37r/ 7r eit 4v / 37r/ 7r e 2 M 



A) 



Li (A) Li(A 2 



A 



A 2 



(3.18) 



A) 



Li (A) Li(A 2 ) 



A 



A 2 



+ 0.84 



c*s(Q 2 

«s(/io) 



2 \ \ 50/81 



A 



A 2 



A 3 



A 4 



(3.19) 



in the case &cleo(x,Q 2 
[D(Q 2 )] res = 



4v / 37r/ 7r ei£ 4v / 37r/ 7r e 2 w 



A 



Li(A 2 ) , L*(A 3 ) L*(A 4 ) 
24 ,\ ; +40 ; n ; -20 



A 2 



A 3 



A 4 



A 



-0.4125 



L/(A) ; A ^ • o,io.-/' U(/i! V" 



A 



A 2 



Li(A) 



A 



«.(Q 2 )\ 



2 x \ 364/405 



a. 



m) 



Li(A) 1on Li(A 



A 



120- 



A 2 



+560 



1112^ + 1008^!) _ 336-^ 



A 4 



A 5 



A 6 



L^A 3 ) 
(3.20) 



also in the case $bms(%, Q 2 ) 
[D(Q 2 )} res = 



Ay/Snf^eJ 4v / 3vr/ 7r e 2 M 



A 



A 



/ Li (A) Li(A 2 ) 



A 



A 2 



+ 0.282 



a s (Q 2 



2 \ \ 50/81 



Li (A) 



A 



Li(A 2 ) Li(A 3 ) Li(A 4 ) , 
24 ,+ ; +40 :„ ; -20 K - ' 1-0.244 



A 2 



A 3 



A 4 



Qa(g2) . 364/405 , {x) {x2) (A3) 



ji 12^1 + 1008 



L*(A 5 ) 



A 



120- 



A 2 



+560- 



336 



L*(A 6 



A 4 A 5 A 6 

where Li (A) is the logarithmic integral for A > 1 defined as the principal value [48] 

dx 



Li (A) = P.V. [ 
Jo 



Inx 



, A = Q 2 /A 2 . 



A 3 

(3.21) 
(3.22) 



Hence, the effective cross section obtained after substitution of the expressions (3.10-3.12) 
into the expression (2.10) is referred as the running coupling effective cross section. We will 
denote the higher-twist cross section obtained using the running coupling constant approach 
by (Err*- 
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IV. CONTRIBUTION OF THE LEADING-TWIST DIAGRAMS 



Regarding the higher-twist corrections to the pion production cross section, a compar- 
ison of our results with leading-twist contributions is crucial. We take two leading-twist 
subprocesses for the pion production: (1) quark- ant iquark annihilation qq — > gj, in which 
the n meson is indirectly emitted from the gluon, g — > it + {it~) and (2) quark-gluon fusion, 
qg — > 57, with subsequent fragmentation of the final quark into a meson, q — > ix + {ix~). The 
corresponding cross sections are obtained in 

d °( qg ^ qi) = - * e >: a ; m (Ui). (4.2) 



dt 3s z \t s 

For the leading-twist contribution, we find 

^ E ^p = ^ / dXldx2dz ( y G ^/hi( x i) G g2/h 2 (x 2 )D^(z)-^^(qq -> # 7 ) + 

Gqi/h 1 {x 1 )G g/h2 (x 2 )D^{z) : ^^{qg -> g 7 )) + * + (4.3) 

where 

Xit „ X 2 M X l t + X 2 U 

s = x±x 2 s, t = , u = , z = . (4.4) 

Z Z X\X 2 S 

Dg(z) = Dg + (z) = Dg ( z) and Dg{z) represents gluon and quark fragmentation functions 
into a meson containing gluon and quark of the same flavor. In the leading-twist subprocess, 
the 7r meson is indirectly emitted from the gluon and quark with the fractional momentum 
z. The 5 function can be expressed in terms of the parton kinematic variables, and the z 
integration can then be done. The final form for the cross section is 

= E— = ^2 dxi dx2 [ G n/hi( x i) G q2/h 2 (x 2 )D n g (z) ■ ——~{qq^gi)+ 

q x lmin J x 2min V 

G qi / hl {x 1 )Gg /h2 (x 2 )D^(z) ■ ~^{ a 9 -»■ 57)^ = 

V- f 1 , f 1 dx 2 ( „ . . „ s D g( z )d<T. 

dx, _ ^G^x^G^xd — ^lM -> 91 ) + 

q u -'- limn v ^2mzn ^ / \ 

XiG qi/hl (x 1 )sx 2 G g / h2 (x 2 ) Dg — ^(qg -> 97)1 • ( 4 - 5 ) 

n at J 
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V. NUMERICAL RESULTS AND DISCUSSION 



In this section, we discuss the numerical results for higher-twist effects with higher-twist 
contributions calculated in the context of the running coupling constant and frozen coupling 
approaches on the dependence of the chosen meson wave functions in the process pp — > 
7r + (or 7r~)7. In the calculations, we use the asymptotic wave function & a sy, the Chernyak- 
Zhitnitsky &cz, the pion wave function (from which two nontrivial Gegenbauer coefficients 
a 2 and 04 have been extracted from the CLEO data on the 7r°7 transition form factor [40]), 
the Braun-Filyanov pion wave functions [41], and the Bakulev-Mikhailov- Stefanis pion wave 
functionf 42]. In Ref.[40], the authors have used the QCD light-cone sum rules approach 
and included into their analysis the NLO perturbative and twist-four corrections. For the 
higher-twist subprocess, we take q\ + q 2 — > (9192) + 7 and we have extracted the following 
four higher-twist subprocesses contributing to pp — > n + (or 7r~)7 cross sections: ud — > 
du — > 7r + 7, ud — > 7T~7, du — * 7r~7 contributing to cross sections. For the dominant leading- 
twist subprocess for the pion production, we take the quark-antiquark annihilation qq — > #7, 
in which the n meson is indirectly emitted from the gluon and quark-gluon fusion, qg — > qj, 
with subsequent fragmentation of the final quark into a meson, q — > it + {it~). As an example 
for the quark distribution function inside the proton, the MRST2003c package [49] has been 
used. The higher twist subprocesses probe the meson wave functions over a large range 
of Q 2 squared momentum transfer, carried by the gluon. Therefore, in the diagram given 
in Fig.l we take Q\ = (xi — l)u, Q\ = —x\t , which we have obtained directly from the 
higher- twist subprocesses diagrams. The same Q 2 has been used as an argument of a s (Q 2 ) 
in the calculation of each diagram. 

The results of our numerical calculations are plotted in Figs. 2-31. First of all, it is very in- 
teresting to compare the resummed higher- twist cross sections with the ones obtained in the 
framework of the frozen coupling approach. In Figs. 2-4 we show the dependence of higher- 
twist cross sections (E^J) calculated in the context of the frozen coupling, (E^ r ) res in the 
context of the running coupling constant approaches and also the ratio R = (E^ r ) res /E^ r )° 
as a function of the pion transverse momentum p T for different pion wave functions at y — 0. 
It is seen that the values of cross sections (E^J) , (E^ r ) res , and R for fixed y and y/s depend 
on the choice of the pion wave function. As seen from Figs. 2-3 in both cases, frozen cou- 
pling and running coupling constant approaches the higher-twist differential cross section is 
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monotically decreasing with an increase in the transverse momentum of the pion. In Figs. 5 
and 6, we shows the dependence of the ratio (Sf + T )°/(S^) and (Sf+ T ) re 7(S^) as a func- 
tion of the pion transverse momentum p T for different pion wave functions. Here (£^ r ) res , 
(S^J) are the higher-twist cross sections calculated in the context of the running coupling 
method and in the framework of the frozen coupling approach and is the leading- 

twist cross section, respectively. As seen from Fig.6, in the region 2 GeV/c < pr < 5 GeV/c 
higher-twist cross section calculated in the context of the running coupling method is sup- 
pressed by about 2 orders of magnitude relative to the leading-twist cross section, but in 
the region 5 GeV/c < px < 30 GeV/c is comparable with the cross section of leading-twist. 
In Figs. 7-10 we show the dependence (Af T )°, (Af T ) res , the ratio r=(Af T ) res / (A% T )° , and 
the ratio (A^ T ) res / (A^ T ) as a function of the pion transverse momentum p? for the pion 
wave functions. Here, (A^ T )° = (Sf + T )° - and (Af r ) rc * = (£f+ T ) res - (£f T ) res . 

As seen from Figs. 7 and 8, the higher-twist differential cross section is decreasing with an 
increase in the transverse momentum of the pion. As is seen from Fig.9, when the transverse 
momentum of the pion is increasing, the ratio r is decreasing. But, as shown in Fig.9, in 
the region 2 GeV / c < pr < 25 GeV / c higher-twist cross section calculated in the context of 
the running coupling method is suppressed by about 3 orders of magnitude relative to the 
higher-twist cross section calculated in the framework of the frozen coupling method. The 
dependence, as shown in Fig. 10, is identically equivalent to Fig.6. 

In Figs. 11-16, we have depicted higher-twist cross sections, ratios (S^ 5 ") , (S^+ r ) res , R = 
(£fJrV(£f7)°, r=(A^rV(A^)°, (Af r )7(A^), and (Af T ) res / (A£ T ) as a function 
of the rapidity y of the pion at = 62.4 GeV and p T = 4.9 GeV/c. At y/s = 62.4 GeV 
and pt = 4.9 GeV jc, the pion rapidity lies in the region —2.52 < y < 2.52. 

As seen from Figs.13-14, in the region (—2.52 < y < —1.92), the ratio for all wave 
functions increase with an increase of the y rapidity of the pion and has a maximum 
approximately at the point y = —1.92. Besides that, the ratio decreases with an in- 
crease in the y rapidity of the pion. As is seen from Figs.13-14, the ratios R and r 
are very sensitive to the choice of the meson wave functions. But, as seen from Fig. 15, 
the ratio (E^) / ) for all wave functions has a minimum approximately at the point 
y = —1.92. In Fig. 16 we show the ratio (A^ T ) res /(A^ T ) as a function of the rapidity y 
of the pion. As seen from Fig. 16, with an increase of the y rapidity of the pion the ra- 
tio increases. It should be noted that the magnitude of the higher-twist cross section for 
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the pion wave function &bms{%, Q 2 ) is very close to the asymptotic wave function$ asy (a;). 
The higher-twist corrections and ratio are very sensitive to the choice of the pion wave 
function. Also, the distinction between R($ asy (x)) with R($cleo( x , Q 2 ))i R(&cz( x ,Q 2 )), 
R(&bf( x ,Q 2 )) and R(§bms( x , Q 2 )) have been calculated. For example, in the case of 
^/s = 62.4 GeV, y = 0, the distinction between R(® asy (x)) with R(§i(x,Q 2 )) (i=CLEO, 
CZ, BF, BMS) as a function of the pion transverse momentum pt is shown in Table [H Thus, 
the distinction between R(Q asy (x)) and Q 2 )){i = CLEO,CZ,BF) is maximum at 

Pt = 20 GeV/c, with R(§ BMS (x)) at pt = 2 GeV/c but the distinction between R(Q asy (x)) 
with R(®i(x, Q 2 )){i = CLEO, CZ, BF) is minimum at p T = 2 GeV/c, with R(®bms( x )) at 
Pt = 20 GeV ' jc and increase with an increase in pt- Such a behavior of R may be explained 
by reducing all moments of the pion model wave functions to those of $> asy (x) for high 
Q 2 . Also, we have calculated the distinction between R(Q asy (x)) with R($ C leo( x , Q 2 )), 
R(§cz( x ,Q 2 )), R(§bf(%,Q 2 )) and R(§bms( x , Q 2 )) as a function of the rapidity y of the 
pion. For example, in the case of y/s = 62AGeV, pt = 4:.9GeV/c the distinction between 
R($asy(x)) with R($i(x,Q 2 )) (i=CLEO, CZ, BF, BMS) as a function of the rapidity y of 
the pion is presented in Table [III 

We have also carried out comparative calculations in the center-of-mass energy -^i = 
200 GeV. The results of our numerical calculations in the center-of-mass energies ^/s = 
200 GeV are plotted in Figs. 17-31. Analysis of our calculations at the center-of-mass energies 
-y/s = 62.4 GeV and -^i = 200 GeV, show that with the increase in beam energy values 
of the cross sections, ratio R = (E^ r ) res /(E^ r )°, and contributions of higher-twist to the 
cross section decrease by about 1-3 order. Therefore the experimental investigation of higher- 
twist effects include renormalon effects conveniently in low energy. On the other hand, the 
higher-twist corrections and ratios R and r are very sensitive to the choice of the pion wave 
function. Also, the distinction between R($ asy (x)) with R($ C leo( x , Q 2 )), R(&cz( x ,Q 2 )), 
R(§bf( x ,Q 2 )) and R(§bms( x , Q 2 )) have been calculated. For example, in the case of 
yfs = 200 GeV, y = 0, the distinction between R($ asy (x)) with R($i(x, Q 2 )) (i=CLEO, CZ, 
BF, BMS) as a function of the pion transverse momentum pt is shown in Table llHl Thus, the 
distinction between R(Q asy (x)) with Q 2 )), (i=CZ, CLEO, BF,) is maximum at pt = 

35 GeV/c, with R(§ B ms( x )) at pt = 10 GeV jc, but the distinction between R(§ asy (x)) 
with R($cz( x , Q 2 )), R($cleo( x , Q 2 )), R($bf( x ,Q 2 )) is minimum at p T = 10 GeV/c, with 
R(§bms{ x )) at Pt = 95 GeV jc and increase with an increase in p T - Also, we have calculated 
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the distinction between R($ asy (x)) with R(§ C leo(x, Q 2 )), R(®cz(x,Q 2 )), R(§bf(x,Q 2 )) 
and R($bms(x, Q 2 ))as a function of the rapidity y of the pion. For example, in the case of 
<fs = 200GeV, p T = Ih.hGeV/c the distinction between R{$ asy (x)) with R($i(x, Q 2 )) as a 
function of the rapidity y of the pion is presented in Table IIVI The calculations show that 
the ratio R($i(x,Q 2 ))/R($ asy (x)), (i=CLEO, CZ, BF, BMS) for all values of the transverse 
momentum px of the pion identically equivalent to ratio r($j(x, Q 2 )) /r(§ asy (x)) . 



VI. CONCLUDING REMARKS 



In this work we have calculated the single meson inclusive production via higher-twist 
mechanism and obtained the expressions for the subprocess qq — > M7 cross section for 
mesons with symmetric wave functions. For calculation of the cross section we have ap- 
plied the running coupling constant method and revealed infrared renormalon poles in the 
cross section expression. Infrared renormalon induced divergences have been regularized by 
means of the principal value prescripton and the resummed expression (the Borel sum) for 
the higher-twist cross section has been found. The higher-twist cross sections were calculated 
in the frozen coupling and running coupling approaches. The resummed higher-twist cross 
section differs from that found using the frozen coupling approach, in some regions, consider- 
ably. Also we demonstrated that higher-twist contributions to single meson production cross 
section in the proton -proton collisions have important phenomenological consequences. We 
have obtained very interesting results. The ratio R for all values of the transverse momen- 
tum pt and of the rapidity of the pion identically equivalent to ratio r. Our investigation 
enables us to conclude that the higher-twist pion production cross section in the proton- 
proton collisions depends on the form of the pion model wave functions and may be used 
for their study. Analysis of our calculations shows that the magnitude of cross sections of 
the leading-twist is larger than the higher-twist cross sections ones calculated in the frozen 
coupling approach in 2-4 order. But, in some regions of transverse momentum of the pion, 
the higher-twist cross section calculated in the context of the running coupling method is 
comparable with the cross sections of leading-twist. Further investigations are needed in 
order to clarify the role of high twist effects in this process. We have demonstrated that 
the resummed result depends on the pion model wave functions used in calculations. The 
proton-proton collisions provide us with a new opportunity to probe a proton's internal 
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structure. In particular, meson production in proton-proton collisions takes into account 
infrared renormalon effects: this opens a window toward new types of parton distributions- 
chiral-odd distributions hi(x,fi 2 ) and h^x, /i 2 ) which can not be measured by the deep 
inelastic lepton-proton scatterings. 
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PT,GeV/c 


R(®cleo(x,Q 2 )) 


R(®cz(x,Q 2 )) 


R(® BF (x,Q 2 )) 


R(^bms{x,Q 2 )) 


U(1?asy(X)) 


H^asy^X)) 


ti[^asy\X)) 


n{<e asy yx)) 


2 


0.557 


0.299 


0.462 


9.813 


6 


1.744 


0.513 


1.5 


2.208 


20 


7.273 


6.065 


6.311 


3.357 



TABLE I: The distinction between R(<$> asy (x)) with R($>i(x,Q 2 )) (i=CLEO, CZ, BF, BMS) at 
cm. energy y/s = 62.4 GeV . 



y 


R(®CLE0(x,Q 2 )) 


R(®cz(x,Q 2 )) 


R($bf(x,Q 2 )) 


R(®bms(x,Q 2 )) 


R(® aS y(x)) 


R{® aS y(x)) 


R(® asy (x)) 


R(®asy(x)) 


-2.52 


11.285 


2.094 


10.162 


3.232 


-1.92 


0.36 


0.279 


0.446 


5.103 


0.78 


0.076 


0.945 


6.213 


2.392 



TABLE II: The distinction between R(<$> asy (x)) with i?($;(x,Q 2 )) (i=CLEO, CZ, BF, BMS) at 
cm. energy y/s = 62.4 GeV and pr = 4.9 GeV/c. 



Pt, GeV jc 


R(®cleo(x,Q 2 )) 

R{$asy(x)) 


R($cz(x,Q 2 )) 

R{$asy(x)) 


R(® BF (x,Q 2 )) 

R($asy(x)) 


R(^bms{x,Q 2 )) 

R($asy(x)) 


10 


1.072 


0.342 


0.625 


2.234 


35 


4.238 


1.115 


4.074 


0.976 


95 


3.011 


0.488 


1,077 


0.561 



TABLE III: The distinction between R(<f> asy (x)) with R(<S>i(x,Q 2 )) (i=CLEO, CZ, BF, BMS) at 
cm. energy y/s = 200 GeV . 
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y 


r(<5>CLEo(x,Q 2 )) 


r(.$cz(x,Q 2 )) 


r(<3> BF (x,Q 2 )) 


R(^bms(x,Q 2 )) 






rCS?asy(X)) 


tiy<i>asy{.X)) 


-2.52 


5.002 


0.823 


3.148 


6.294 


-1.92 


1.538 


0.285 


0.447 


1.089 


0.78 


0.504 


0.861 


4.351 


4.149 



TABLE IV: The distinction between R(<$> asy (x)) with R(<S>i(x,Q 2 )) (i=CLEO, CZ, BF, BMS) at 
cm. energy y/s = 200 GeV and px = 15.5 GeV/c. 
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FIG. 1: Feynman diagrams for the higher-twist subprocess, q\q<2. — ► vr + (or tt ) r y. 
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FIG. 2: Higher-twist tt + production cross section (T, HT )° as a function of the pt transverse mo- 
mentum of the pion at the cm. energy \fs = 62.4 GeV. 




p T , GeV/c 



FIG. 3: Higher-twist 7r + production cross section (Y l HT ) res as a function of the px transverse 
momentum of the pion at the cm. energy y/s = 62.4 GeV. 
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FIG. 4: Ratio R = (£^ r ) re V(^^+ r )°> where higher-twist contribution are calculated for the pion 
rapidity y = at the c.m.energy \/s = 62.4 GeV as a function of the pion transverse momentum, 

PT- 




FIG. 5: Ratio ) / (S^+ ), as a function of the pt transverse momentum of the pion at the cm. 
energy y/s = 62.4 GeV. 
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FIG. 6: Ratio ) res /(S^+ ), as a function of the px transverse momentum of the pion at the 
cm. energy ^fs = 62.4 GeV. 




FIG. 7: The difference of the higher-twist cross section, (A^ T )° = (SfJ) - as a function 

of the pion transverse momentum, px, at the cm. energy y/s = 62.4 GeV. 
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FIG. 8: The difference of the higher-twist cross section, (A* T ) res = (Sf+ T ) res - {^"-) res , as a 
function of the pion transverse momentum, pt, at the c.m.energy \fs = 62.4 GeV. 
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FIG. 9: Ratio r = (A^ ) res /(A^ T )°, where higher- twist contributions are calculated for the pion 
rapidity y = at the cm. energy y/s = 62.4 GeV, as a function of the pion transverse momentum, 
Pt- 
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FIG. 10: Ratio (A^f ) res /(A^ ), where higher-twist contributions are calculated for the pion 
rapidity y = at the cm. energy y^i = 62.4 GeV, as a function of the pion transverse momentum, 
PT- 
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FIG. 11: Higher- twist 7r + production cross section (T,^+)°, as a function of the y rapidity of the 
pion at the transverse momentum of the pion px = 4.9 GeV/c, at the cm. energy ^/s = 62.4 GeV. 
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FIG. 12: Higher-twist 7r + production cross section (S^ r ) res , as a function of the y rapidity of the 
pion at the transverse momentum of the pion = 4.9 GeV/c, at the cm. energy ^/s = 62.4 GeV. 
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FIG. 13: Ratio R = (£^ r ) res /(£^ T )°, as a function of the y rapidity of the pion at the transverse 
momentum of the pion pt = 4.9 GeV/c, at the cm. energy y/s = 62.4 GeV. 
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FIG. 14: Ratio r = (A^ T ) res / (A^ T )° , as a function of the y rapidity of the pion at the transverse 
momentum of the pion pt = 4.9 GeV/c, at the cm. energy y/s = 62.4 GeV. 




FIG. 15: Ratio (S^ r )°/(S^J), as a function of the y rapidity of the pion at the transverse mo- 
mentum of the pion = 4.9 GeV/c, at the cm. energy yfs = 62.4 GeV . 
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FIG. 16: Ratio (S^ r ) res /(S^), as a function of the y rapidity of the pion at the transverse 
momentum of the pion pr = 4.9 GeV/c, at the cm. energy y/s = 62 A GeV. 




FIG. 17: Higher- twist tt + production cross section )° as a function of the pr transverse 

momentum of the pion at the cm. energy \/s = 200 GeV. 
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FIG. 18: Higher-twist ir + production cross section ) res as a function of the pt transverse 

momentum of the pion at the cm. energy \/s = 200 GeV. 




FIG. 19: Ratio R = (S^ r ) res /(S^ r )°, where higher-twist contribution are calculated for the pion 
rapidity y = at the cm. energy y/s = 200 GeV as a function of the pion transverse momentum, 
Pt- 
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FIG. 20: Ratio (S^)°/(E^+), as a function of the px transverse momentum of the pion at the 
cm. energy ^fs = 200 GeV . 
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FIG. 21: Ratio (S^ r ) res /(E^+ ), as a function of the pt transverse momentum of the pion at the 
cm. energy y/s = 200 GeV. 
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FIG. 22: The difference of the higher-twist cross section, (A^ T )° = (£f J)° - (X^) , as a function 
of the pion transverse momentum, px, at the cm. energy y/s = 200 GeV. 
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FIG. 23: The difference of the higher-twist cross section, (A^ T ) res = (£f+ T ) res - (Sf„ T ) res , as a 
function of the pion transverse momentum, px, at the c.m.energy y/s = 200 GeV. 
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FIG. 24: Ratio r = (A^ T ) res /(A^ )°, where higher- twist contributions are calculated for the pion 
rapidity y = at the cm. energy y/s = 200 GeV, as a function of the pion transverse momentum, 
PT- 




FIG. 25: Ratio (A^ T ) res /(A^ T ), where higher-twist contributions are calculated for the pion 

rapidity y = at the cm. energy ^/s = 200 GeV, as a function of the pion transverse momentum, 
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FIG. 26: Higher-twist ir + production cross section (T,^) , as a function of the y rapidity of the 
pion at the transverse momentum of the pion px = 15.5 GeV/c, at the cm. energy y/s = 200 GeV. 
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FIG. 27: Higher-twist 7r + production cross section ) res , as a function of the y rapidity of the 
pion at the transverse momentum of the pion pt = 15.5 GeV/c, at the cm. energy ^fs = 200 GeV. 
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FIG. 28: Ratio R = (S^ r ) res /(S^ T )°, as a function of the y rapidity of the pion at the transverse 
momentum of the pion pt = 15.5 GeV/c, at the cm. energy y/s = 200 GeV. 




FIG. 29: Ratio r = (A^ T ) res / (A^ T )° , as a function of the y rapidity of the pion at the transverse 
momentum of the pion pt = 15.5 GeV/c, at the cm. energy y/s = 200 GeV. 
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FIG. 30: Ratio (S^ r )°/(S^+ ), as a function of the y rapidity of the pion at the transverse mo- 
mentum of the pion p? = 15.5 GeV/c, at the cm. energy yfs = 200 GeV. 
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FIG. 31: Ratio ) res /(S^+ ), as a function of the y rapidity of the pion at the transverse 

momentum of the pion pt = 15.5 GeV/c, at the cm. energy y/s = 200 GeV. 
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